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i.e.,  (a)  contact-pressure  history,  (h)  reciprocation  frequency,  and  (c)  reciprocation  amplitude.  Examination 
of  the  results  obtained  clearly  revealed  the  presence  of  three  zones  within  the  weld,  i.e.,  (a)  Contact-interface 
region,  (b)  Thermo-mechanically  affected  zone,  and  (c)  heat-affected  zone.  While  there  are  no  publicly 
available  reports  related  to  Carpenter  Custom  465  LFW  behavior,  to  allow  an  experiment/computation 
comparison,  these  findings  are  consistent  with  the  results  of  our  ongoing  companion  experimental 
investigation. 
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1.  Introduction 

Linear  friction  welding  (LFW)  falls  into  the  category  of 
solid-state  frictional-sliding-based  joining  processes.  As  sche¬ 
matically  shown  in  Fig.  1,  the  LFW  process  involves  (relative) 
reciprocating  linear  motion  of  the  two  workpieces  to  be  joined. 
As  a  result  of  frictional  sliding,  heat  is  dissipated  at  the  faying/ 
contacting  flat  surfaces,  promoting  material  softening  in  the 
surrounding  region,  plastic  deformation,  material  expulsion, 
and  ultimate  adhesion/bonding. 

Detailed  past  investigations  of  LFW  (Ref  1-3)  established  the 
presence  of  five  distinct  phases  of  this  joining  process:  (a)  initial 
phase,  (b)  transition  phase,  (c)  equilibrium  phase,  (d)  deceleration/ 
forging  phase,  and  (e)  stand  phase.  A  simple  schematic  of  the  first 
four  phases  is  depicted  in  Fig.  2.  Since  the  last  phase  does  not 
involve  workpiece  length-scale  structural/kinematic  effects,  it  is  not 
shown.  Considering  the  fact  that  the  defining  features  of  the  five 
LFW  phases  can  be  found  in  Ref  1  -3,  they  will  not  be  repeated  here. 
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The  main  LFW  process  parameters  can  be  defined  as 
follows:  (a)  alloy-grade  and  thermal-mechanical  treatment  of 
the  workpiece  materials  to  be  joined,  (b)  frequency  of 
reciprocating  motion,  (c)  amplitude  of  reciprocating  motion, 
(d)  temporal  variation  of  contact/upsetting/forging  pressure,  (e) 
total  duration  of  LFW  process,  and  (f)  the  geometry  of  the 
contacting  surfaces.  It  is  generally  accepted  that  selection  of 
these  parameters  (or  more  precisely,  their  combination)  may 
have  a  profound  effect  on  the  structural  integrity  and  quality  of 
the  resulting  weld. 

A  detailed  account  of  the  main  advantages  of  LFW  relative 
to  the  other  friction-based  welding  processes,  solid-state  joining 
processes,  and  fusion-based  welding  processes  can  be  found  in 
Ref  4.  The  same  reference  also  provides  a  detailed  description 
of  the  key  limitations,  as  well  as  the  main  areas  of  application, 
of  this  joining  process. 

Examination  of  the  open-domain  literature  carried  out  as 
part  of  the  present  work  identified  a  number  of  investigations 
demonstrating  successful  application  of  the  LFW  process  to 
join  various  metal-alloy  grades  such  as  (a)  titanium  alloys  [e.g., 
(Ref  5-8)],  (b)  superalloys  (Ref  9),  (c)  steels  (Ref  10,  11),  and 
(d)  intermetallics  (Ref  12).  Cumulatively,  these  investigations 
established  that  the  temperature  in  the  contact  region  plays  a 
dominant  role  in  the  joining  process,  affecting  the  phenomena 
such  as  welding  efficiency,  spatial  distribution,  and  temporal 
evolution  of  material  deformation  and  micro  structure  fields,  as 
well  as  quality/structural  integrity  of  the  weld. 

The  work  presented  in  this  manuscript  deals  with  LFW  of 
Carpenter  Custom  465  precipitation-hardened  martensitic  stain¬ 
less  steel.  This  material  is  being  considered  for  use  in  aircraft- 
structure  applications,  and  LFW  is  being  considered  as  the  main 
joining-process  candidate.  However,  no  reports  pertaining  to 
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Reciprocating  Workpiece 


Fig.  1  A  schematic  of  a  single  cycle  of  the  Linear  Friction  Welding  (LFW)  process 


Fig.  2  First  four  phases  of  the  LFW  process:  (a)  initial  phase,  (b)  transition  phase,  (c)  equilibrium  phase,  and  (d)  deceleration/forging  phase 


LFW  of  this  material  could  be  found  in  the  open  literature.  The 
present  work  tries  to  fill  this  void  through  the  development, 
parameterization,  and  validation  of  an  LFW  process  model. 


The  material  closest  to  Carpenter  Custom  465  for  which  an 
experimental  LFW  investigation  has  been  reported  is  AISI 
3 16L,  a  low-carbon  austenitic  stainless  steel  (Ref  13).  The  work 
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presented  in  Ref  13  established  three  important  points:  (a) 
through  proper  selection  of  the  LFW  process  parameters, 
sound/flaw-free  welds  can  be  produced;  (b)  again  through 
proper  selection  of  the  process  parameters,  the  development  of 
undesirable  micro  structures  within  different  weld  zones  can  be 
avoided.  For  example,  by  selecting  the  process  parameters  to 
maximize  the  bum-off  rate  (i.e.,  the  rate  at  which  plasticized 
material  is  expelled  from  the  contact  region),  formation  of  the 
undesirable  5-ferrite  can  be  prevented  in  the  weld  region;  and 
(c)  determination  of  the  optimum  process  parameters  using 
purely  experimental  means  relies  on  the  post-mortem  (typically 
destmctive)  examination  of  the  LFW  joints. 

As  established  above,  a  purely  experimental  approach  to  the 
LFW  process  suffers  from  a  number  of  shortcomings  related  to 
the  real-time  determination  of  the  (optimal)  process  parameters 
and  the  weld  conditions  (e.g.,  temperature  distribution  over  the 
contact  surface  and  in  the  direction  normal  to  it),  micro  structure 
distribution  throughout  the  weld  region,  etc.  These  shortcom¬ 
ings  can  be  potentially  overcome  through  the  use  of  high- 
fidelity  LFW  process  models.  An  overview  of  the  public- 
domain  literature  revealed  the  following  LFW  modeling  and 
simulation  studies  which  appear  the  most  noteworthy:  (a)  a 
two-dimensional  fully  coupled  thermo-mechanical  finite-ele¬ 
ment  model  of  the  LFW  process  has  been  developed  by 
D’Alvise  et  al.  (Ref  14),  who  employed  an  innovative  adaptive 
contact  and  re-meshing  algorithm  to  simulate  flash  formation. 
However,  due  to  the  two-dimensional  character  of  the  process 
model,  full  validation  of  the  model  against  the  available 
experimental  data  was  associated  with  a  number  of  challenges; 
(b)  a  three-dimensional  fully  coupled  thermo-mechanical  finite- 
element  model  for  the  LFW  process  has  been  developed  by 
Sorina-Muller  et  al.  (Ref  15),  who  applied  it  to  Ti-6Al-2Sn- 
4Cr-6Mo  alloy  weldments.  While  the  computed  results  and  the 
corresponding  experimental  data  related  to  the  temporal 
evolution  of  the  axial  shortening  and  the  spatial  distribution 
of  microhardness  in  the  as-welded  material  state  are  claimed  to 
be  in  very  good  agreement,  the  constitutive  model  for  the 
subject  material  was  fairly  simple  and  did  not  properly  reflect 
many  aspects  of  the  material  mechanical  response  under 
dynamic,  elevated-temperature,  and  high-contact-pressure  load¬ 
ing  conditions;  and  (c)  the  model  of  Sorina-Muller  et  al.  (Ref 
15)  has  been  extended  by  Grujicic  et  al.  (Ref  16)  in  order  to 
remove  the  aforementioned  material-model  shortcomings  and 
provide  a  deeper  insight  into  the  development  of  material 
microstructure  (and  the  associated  mechanical  properties) 
within  different  weld  regions  and  their  dependence  on  the 
LFW  process  parameters.  Specifically,  since  the  HAZ  of  the 
weld  has  been  found  to  possess  somewhat  inferior  mechanical 
properties  relative  to  the  TMAZ  and  the  base  metal,  particular 
attention  has  been  given  to  the  problem  of  establishing 
processing/microstructure/property  relations  in  this  weld  re¬ 
gion.  The  model  predictions  and  their  experimental  counter¬ 
parts  pertaining  to  the  overall  structural  performance  of  the 
LFW  joint  are  found  to  be  in  reasonably  good  agreement. 

The  main  objective  of  the  present  work  is  to  develop, 
parameterize,  and  validate  an  LFW  process  model  for  Carpen¬ 
ter  Custom  465.  Since  a  thermo-mechanical  constitutive  model 
for  this  material  is  not  presently  available  in  the  open  literature, 
a  major  effort  in  this  work  will  be  directed  toward  the 
development,  parameterization,  and  validation  of  the  appropri¬ 
ate  material  model.  In  addition,  since  contact  interactions  of  the 
workpieces  to  be  joined  play  a  key  role  in  the  LFW-joint 
formation  and  the  accompanying  changes  in  the  material 


micro  structure,  considerable  effort  will  also  be  directed  toward 
improving  the  contact  models  and  their  parameterization. 


2.  Carpenter  Custom  465 

As  mentioned  earlier,  the  material  analyzed  in  the  present 
work  is  the  premium  (i.e.,  double  vacuum)-melted  Carpenter 
Custom  465  precipitation-hardened  martensitic  stainless  (i.e., 
oxidation-  and  corrosion-resistant)  steel  (Ref  17).  In  this  steel, 
improved  oxidation  resistance  is  achieved  through  the  addition 
of  chromium  (minimum  ca.  12  wt.%)  which  ensures  formation 
of  a  continuous  Cr203-based  oxide  protective  surface  film  (with 
low-oxygen  permeability  and  high  adhesion  strength).  En¬ 
hanced  corrosion  resistance,  on  the  other  hand,  is  achieved  by 
ensuring,  through  the  chemical  composition  and  heat  treatment, 
that  the  matrix  contains  only  one  crystalline  phase  (i.e., 
martensite,  in  the  present  case).  In  (precipitation-hardened) 
Carpenter  Custom  465,  the  requirement  for  increased  fracture 
toughness  and  for  the  enhanced  retention  of  mechanical 
properties  at  high  temperatures  demands  the  introduction  of  a 
secondary/precipitate  phase  (using  an  aging  heat  treatment).  To 
ensure  high  resistance  toward  “dissimilar-metal”  corrosion 
(i.e.,  corrosion  between  the  matrix  and  the  precipitate  phases), 
the  two  phases  are  designed  to  have  a  small  electronegativity 
difference.  This  is  achieved  by  ensuring  that  the  precipitate 
phase  is  of  an  intermetallic  character  (Ref  18). 

2.1  Heat  Treatment(s) 

Carpenter  Custom  465  with  the  temper  designation  H900  or 
HI 000  is  typically  used  in  aircraft  frames.  It  is,  hence,  the  LFW 
of  these  Carpenter  Custom  465  tempers  that  is  analyzed  in  the 
present  work.  The  H900/1000  heat  treatment  involves  the 
following  three  steps: 

(a)  Solutionizing  heat  treatmen — Within  this  step,  the  mate¬ 
rial  is  first  subjected  to  high  temperatures  (e.g.,  1255  K/ 
1  h),  to  convert  its  as-received  microstructure  into  aus¬ 
tenite,  and  then  quenched  to  room  temperature  (to  pre¬ 
vent  diffusion-type  decomposition  of  austenite).  Since 
the  room  temperature  is  between  the  temperatures  for 
the  start  and  finish  of  the  martensitic  transformation  (the 
so-called  Ms  and  Mf  temperatures,  respectively),  the  as- 
quenched  material  consists  of  untransformed  austenite 
and  martensite; 

(b)  Cold-Treating — In  this  step,  the  as-quenched  material  is, 
within  a  time  period  of  less  than  24  h  following  the  so¬ 
lutionizing  treatment,  subjected  to  a  refrigeration  treat¬ 
ment  at  200  K  for  about  eight  hours.  Since  Mf  is  higher 
than  this  temperature,  austenite-to-martensite  conversion 
continues  until  (near)  completion.  Subsequently,  the 
material  is  allowed  to  warm  to  room  temperature.  The 
resulting  fully  martensitic  structure  contains  numerous 
lath  interfaces  and  intra-lath  dislocations,  and  so  pro¬ 
vides  an  abundance  of  potential  nucleation  sites  for  the 
precipitate  phase(s).  This,  in  turn,  ensures  that  during 
the  subsequent  aging  treatment,  the  resulting  microstruc¬ 
ture  will  contain  a  high  number  density  of  very  fine 
(material-property-enhancing)  precipitates;  and 

(c)  Aging  Heat  Treatment — In  this  step,  the  annealed, 
quenched,  and  refrigeration-treated  material  is  aged  in  a 
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temperature  range  between  755  and  894  K  for  four  to 
eight  hours.  This  is  followed  by  a  water  or  oil  quench 
of  the  material.  The  resulting  as-aged  material  condition 
is  designated  as  Hxxx,  where  xxx  is  replaced  with  the 
aging  temperature  expressed  in  °F  (e.g.,  HI 000).  A 
trade-off  in  material  strength,  fracture  toughness,  and 
stress-corrosion-cracking  resistance  in  aged  Carpenter 
Custom  465  (Ref  19),  caused  by  differences  in  the  aging 
temperature,  are  demonstrated  in  Fig.  3(a-c).  It  is  seen 
that  as  the  aging  temperature  increases,  material  strength 
decreases,  Fig.  3(a),  while  its  fracture  toughness, 
Fig.  3(b),  and  stress-corrosion-cracking  resistance, 
Fig.  3(c),  increase.  These  results  suggest  that  the  expo¬ 
sure  of  the  material  within  the  weld  to  high  temperatures 
can  result  in  over-aging  and,  thus,  in  the  (generally 
undesirable)  loss  of  strength. 


Age-Hardening  Temperature,  deg.  F 
900  1000  1100 


900  1000  1100 


22  Precipitation  Reaction(s)  During  Aging 

In  the  course  of  the  aging  treatment  of  Carpenter  Custom 
465,  coherent  hexagonal- structure  needle-shaped  (material- 
strength-controlling)  co-phase  precipitates  are  first  formed.  As 
the  aging  temperature/duration  is  increased,  and  the  material 
begins  to  over-age,  these  precipitates  gradually  become 
replaced  by  incoherent  orthorhombic-structure  plate-like 
Ni3(Ti,Mo)  precipitates.  These  micro  structural  changes  are 
accompanied  by  a  loss  of  material  strength  and,  initially,  an 
improvement  in  its  fracture  toughness.  At  sufficiently  high 
aging  temperatures,  martensite-to-austenite  reversion  may  take 
place,  giving  rise  to  a  severe  loss  of  material  strength  (a  highly 
undesirable  effect)  and  an  increase  in  ductility.  Changes  in  the 
micro  structure  of  the  material  within  the  weld  in  Carpenter 
Custom  465  are  currently  being  investigated  in  our  ongoing 
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Fig.  3  The  effect  of  aging  temperature  on  the  material  (a)  yield  strength,  (b)  fracture  toughness,  and  (c)  stress-corrosion-cracking  (SCC)  resis¬ 
tance  at  room  temperature  in  Carpenter  Custom  465 
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work.  In  the  present  manuscript,  the  focus  will  be  placed  on 
determining  the  thermo-mechanical  conditions  within  the  weld 
region  which  cause  these  changes 


3.  LFW  Process  Modeling 

3.1  Problem  Definition 

As  mentioned  earlier,  the  main  objective  of  the  present  work 
is  to  develop  a  process  model  which  can  be  used  to  investigate 
the  effect  of  the  LFW  process  parameters  in  Carpenter  Custom 
465  on  the  spatial  distribution  and  temporal  evolution  of 
temperature  and  material  flow  in  the  weld  region,  as  well  as  the 
temporal  evolution  of  axial  shortening.  In  our  ongoing  work, 
this  model  is  being  extended  in  order  to  enable  prediction  of  the 
material  micro  structural  changes  in  different  LFW  weld  zones. 

3.2  Formulation  of  the  Process  Model 

Modeling  of  the  LFW  process  carried  out  in  the  present 
work  employed  a  fully  coupled  thermo-mechanical  finite- 
element  procedure,  which  is  an  adaptation  and  extension  of  the 
one  developed  for  and  applied  to  the  friction  stir  welding 
(FSW)  process  in  our  prior  work  (Ref  20-28).  The  approach 
used  requires  specification  of  the  following:  (a)  computational 
domain;  (b)  computational-analysis  type;  (c)  initial  conditions; 
(d)  boundary  conditions;  (e)  contact  interactions;  (f)  heat- 
generation  and  partitioning;  (g)  computational  algorithm;  (h) 
material  thermo-mechanical  models;  and  (i)  computational 
accuracy,  stability,  and  cost.  These  aspects  are  briefly  over¬ 
viewed  in  the  remainder  of  this  sub-section. 

3.3  Computational  Domain 

The  computational  domain  used  consists  of  two  (initially) 
parallelepiped-shaped  workpieces  with  edge  dimensions  Lx  by 
Ly  by  Lz  =  20  mm  by  10  by  30  mm.  The  two  workpieces  are 
stacked  in  the  z-direction  so  that  the  initial  cross-sectional  area 
of  the  two  faying  surfaces  is  Lx  by  Ly  (=200  mm2).  In 
accordance  with  the  common  LFW  practice,  reciprocating 
motion  was  assumed  to  take  place  along  the  shorter  of  the  two 
contact-interface  directions  (y-direction,  in  the  present  case). 
This  is  done  in  order  to  facilitate  flash  formation  (i.e.,  expulsion 
of  the  plasticized  material)  and  to  promote  the  formation  of  a 
good  quality  weld. 

Each  of  the  workpieces  is  meshed  using  first-order  eight-node 
hexahedron  (cubic  or  nearly  cubic)  thermo-mechanically  cou¬ 
pled,  reduced-integration  finite  elements.  Due  to  the  transient 
nature  of  the  temperature  field  and  high-temperature  sensitivity 
of  the  workpiece  material  properties,  mesh  size  (in  the  direction 
normal  to  the  faying  surfaces)  is  chosen  in  accordance  with  the 
steepest  temperature  gradient  (and  the  stable  time  increment 
desired).  Due  to  the  highly  localized  nature  of  the  HAZ  and 
TMAZ,  a  fine  mesh  had  to  be  used  only  in  the  workpiece  region 
near  the  contact  interface.  Typically,  each  workpiece  initially 
contained  ca.  20,000  elements.  An  example  of  a  typical  finite- 
element  mesh  used  is  depicted  in  Fig.  4  in  which,  for  clarity,  only 
the  lower  (stationary)  workpiece  is  shown. 

3.4  Computational-Analysis  Type 

The  LFW  process  is  analyzed  computationally  using  a  fully 
coupled  thermo-mechanical  finite-element  algorithm,  within 


which  heat  dissipation  associated  with  workpiece/workpiece 
interfacial  friction-sliding  and  plastic  deformation  is  treated  as  a 
heat  source  in  the  governing  heat  conduction  equation,  while 
the  effect  of  temperature  on  the  mechanical  response  of  the 
workpiece  material  is  taken  into  account  through  the  use  of 
temperature-dependent  material-model  parameters. 

3.5  Initial  Conditions 

The  analysis  is  carried  out  by  prescribing,  from  the  onset, 
zero  values  to  the  stress,  and  material-particle  velocities  and  the 
ambient  value  to  the  temperature  in  both  workpieces. 

3.6  Boundary  Conditions 

While  LFW  requires  only  the  relative  reciprocating  motion 
of  the  contacting  upper  and  lower  workpieces,  for  convenience, 
the  lower  workpiece  is  kept  stationary  in  the  present  work  while 
the  reciprocating  motion  is  assigned  to  the  upper  workpiece. 
This  is  accomplished  by  assigning  zero-displacement  boundary 
conditions  (in  all  three  directions)  to  the  bottom  face  of  the 
lower  workpiece.  As  far  as  the  top  face  of  the  upper  workpiece 
is  concerned,  it  is  constrained  in  the  x-direction,  subjected  to  a 
time-dependent  pressure  in  the  z-direction  and  assigned  a 
sinusoidal/reciprocating  displacement  in  the  y-direction,  Uy,  as 

uy  =  uo  sin  (2nft),  (Eq  1) 

where  /  is  the  (ordinary)  frequency  of  reciprocating  motion 
(Hz)  and  u0  is  the  reciprocation  amplitude. 

As  far  as  the  thermal  boundary  conditions  are  concerned,  it 
is  assumed  that  the  temperature  at  the  outer  surfaces  of  the 
workpieces  is  controlled  by  forced  convection  conditions 
(induced  by  the  relative  reciprocating  motion  of  the  workpiec¬ 
es)  and  by  radiation  effects.  Typical  values  for  the  forced 
convection  coefficient  of  100  W/m2K  and  for  the  emissivity  of 
0.5  were  chosen.  The  same  boundary  conditions  were  applied 
to  the  portions  of  the  faying  surfaces  which  were  exposed 
(temporarily)  as  a  result  of  the  reciprocating  motion  of  the 
workpieces. 

3.7  Contact  Interactions 

The  workpiece-workpiece  normal  interactions  are  analyzed 
using  a  penalty-contact  algorithm.  Within  this  algorithm, 
(normal)  penetration  of  the  contacting  surfaces  is  resisted  by 
a  set  of  linear  springs  which  produce  a  contact  pressure  that  is 
proportional  to  the  depth  of  penetration.  Typically,  maximum 
default  values,  which  still  ensure  computational  stability,  are 
assigned  to  the  (penalty)  spring  constants.  Force  equilibrium  in 
a  direction  collinear  with  the  contact-interface  normal  then 
causes  the  penetration  to  acquire  an  equilibrium  (contact- 
pressure-dependent)  value.  It  should  be  noted  that  no  contact 
pressures  are  developed  unless  (and  until)  the  nodes  on  the 
slave  surface  contact/penetrate  the  master  surface.  On  the  other 
hand,  the  magnitude  of  the  contact  pressure  that  can  be 
developed  is  unlimited.  As  far  as  the  tangential  workpiece- 
workpiece  interactions  (responsible  for  transmission  of  the 
shear  stresses  across  the  contact  interface)  are  concerned,  they 
are  modeled  using  a  modified  Coulomb  friction  law.  Within  this 
law,  the  maximum  value  of  the  shear  stresses  that  can  be 
transmitted  (before  the  contacting  surfaces  begin  to  slide)  is 
defined  by  a  product  of  the  contact  pressure  and  a  static  (before 
sliding)  or  a  kinetic  (during  sliding)  friction  coefficient.  In 
addition,  to  account  for  the  potential  occurrence  of  a  sticking 
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Fig.  4  Typical  finite-element  mesh  used  in  the  present  work.  For  clarity,  only  the  lower  workpiece  is  shown 


condition  (sliding  occurs  by  shear  fracture  of  the  softer 
workpiece  material,  in  the  case  of  LFW  of  dissimilar  materials, 
rather  than  by  a  relative  motion  at  the  contact  interface),  a 
maximum  value  of  shear  stress  (equal  to  the  temperature- 
dependent  shear  strength  of  the  softer  material)  that  can  be 
transmitted  at  any  level  of  the  contact  pressure,  is  also 
specified. 

As  far  as  the  friction  coefficient  is  concerned,  it  is  generally 
assumed  that  in  the  case  of  friction-based  welding  processes, 
this  contact  parameter  is  controlled  by  the  formation  and 
shearing  of  micro  welds  (i.e.,  micron-sized  regions  at  which 
contacting  surface  asperities  are  bonded).  Furthermore,  it  is 
recognized  that  the  friction  coefficient  is  a  function  of  a  number 
of  factors  such  as  the  contact  interface  (mean)  temperature,  slip 
speed,  contact  pressure,  faying-surfaces’  roughness/topology, 
etc.  Using  the  experimental  and  numerical  results  pertaining  to 
the  effect  of  various  contact-interface  conditions  on  the  friction 
coefficient,  reported  in  Ref  29,  prototypical  functional  relations 
between  steel/steel  friction  coefficient  and  temperature,  slip 
velocity,  and  contact  pressure  are  depicted  in  Fig.  5(a-c), 
respectively.  The  functional  relations  depicted  in  these  figures 
reflect  the  effect  of  the  three  contact-surface  variables  on  the 
competition  between  microweld  formation  and  shear-fracture 
processes.  Specifically,  increased  temperature  and  sliding  rate 
lead  to  material  softening,  which  both  increases  the  tendency 
for  micro  weld  formation  and  for  shear  fracture.  As  far  as  the 
effect  of  contact  pressure  is  concerned,  it  can  promote  the 
development  of  a  local  texture  within  the  microwelds  and,  in 
turn,  a  low  value  of  the  directionally  dependent  shear  strength. 

3.8  Heat-Generation  and  Partitioning 

As  mentioned  earlier,  both  frictional  sliding  and  plastic 
deformation  act  as  heat  sources  during  the  LFW  process.  The 
heat  generated  by  these  two  phenomena  in  the  workpiece- 
contact  region  is  subsequently  partitioned  between  the  two 
workpieces. 

To  account  for  the  heat  generated  due  to  frictional  sliding,  its 
rate  is  assumed  to  scale  with  the  product  of  local  interfacial 
shear  stress  and  the  sliding  rate. 

As  far  as  work  due  to  plastic  deformation  is  concerned,  95% 
of  this  work  was  assumed  to  be  dissipated  in  the  form  of  heat, 


while  the  remaining  5%  is  stored  in  the  form  of  crystalline 
defects. 

Partitioning  of  the  heat  generated  during  LFW  between  the 
two  workpieces  is  computed  using  the  procedure  described  in 
Ref  30.  According  to  this  procedure,  a  ratio  of  the  fractions  of 
frictional-sliding-generated  heat  allotted  to  the  upper  and  lower 
workpieces  is  given  by 
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(Eq  2) 


where  X  is  thermal  conductivity,  p  is  the  mass  density,  and  cp 
is  the  specific  heat  of  the  workpiece  material.  It  should  be 
noted  that  the  (q upper  +  Slower)  sum  must  be  equal  to  the  total 
heat  generated  during  frictional  sliding.  Examination  of  Eq 
(2)  reveals  that  for  identical  materials  of  the  two  workpieces, 
the  heat  ratio  is  equal  to  unity  as  expected.  Since  Eq  (2)  is 
frequently  used  in  computational  analyses  of  various  fric- 
tional-sliding-based  welding  processes,  e.g.,  LFW,  FSW,  etc., 
without  providing/explaining  the  basic  foundation  of  this 
equation,  such  an  explanation  is  provided  in  the  remainder  of 
this  paragraph.  To  rationalize  Eq  (2),  the  following  observa¬ 
tions  should  be  made:  (a)  at  any  time  t ,  the  distance  traveled 
by  the  heat  front  from  the  contact  interface  within  each  of  the 
workpieces  scales  with  y'a,  where  a  is  the  respective  material 
thermal  diffusivity;  (b)  for  a  given  instantaneous  value  of  the 
interfacial  temperature,  temperature  profile  integrated  over  the 
distance  traveled  by  the  heat  wave  from  the  interface  at  a  gi¬ 
ven  time  also  scales  with  y'a;  and  (c)  heat  absorbed  by  the 
workpiece  material  in  the  time  period  0  to  t  scales  with  the 
product  of  the  material  mass  density  p,  its  mass-based  spe¬ 
cific  heat  cp,  and  the  aforementioned  integral,  i.e.,  with 
pcp^/a.  Furthermore,  since  a  is  related  to  thermal  conductiv¬ 
ity  X  as  a  =  X/Xpcp.pcp,  the  heat  absorbed  scales  with 
\/^P  cp ■ 

In  addition  to  the  heat  partitioning  between  the  contacting 
workpieces,  the  problem  of  heat  transfer  between  the  work- 
pieces  had  to  be  addressed.  In  the  present  work,  this  heat 
transfer  was  assumed  to  be  purely  conductive.  In  other  words, 
while  in  the  initial  phase  of  LFW  the  two  workpieces  contact 
only  through  surface  asperities  and  air-filled  pockets  exist  in  the 
contact  region,  the  contribution  of  convection  and  radiation  to 
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Fig.  5  Variation  of  friction  coefficient  with  (a)  temperature,  (b)  slip  rate,  and  (c)  contact  pressure  in  prototypical  steel/steel  sliding  contact  pairs 
(Ref  29) 


the  heat  exchange  across  these  pockets  was  deemed  negligible. 
This  was  done  for  the  following  reasons:  (a)  in  the  initial  LFW 
stage,  the  temperature  in  the  contact  region  is  too  low  for 
radiation  to  play  any  significant  role;  and  (b)  the  interfacial 
pockets  are  too  small  and  isolated  for  a  convection  current  to  be 
formed.  As  far  as  the  conductive  heat  transfer  between  the 
contacting  surfaces  is  concerned,  it  is  defined  as  the  product  of 
the  temperature  difference  between  the  contacting  surfaces  and 
the  interfacial  thermal  conductance.  It  is  well  established  that 
the  interfacial  thermal  conductance  is  a  function  of  the  surface 
topology,  contact  pressure,  and  material  thermo-mechanical 
properties  (Ref  31).  An  example  of  typical  ( T  <  500  K,  in  air) 
results  pertaining  to  the  effect  of  contact  pressure  and  mean 
surface  roughness  on  the  stainless-steel/stainless-steel  interfa¬ 
cial  thermal  conductance  is  depicted  in  Fig.  6  (Ref  32).  Under 
high-contact  pressures  and  high  interfacial  temperatures 
encountered  in  the  course  of  LFW,  plastic  deformation  of  the 
surface  asperities  gives  rise  to  high  values  of  the  interfacial 
thermal  conductance.  Consequently,  temperature  jumps  across 


the  contact  interface  are  expected  to  be  very  small,  and  are 
ignored  in  the  present  work. 

3.9  Computational  Algorithm 

As  established  earlier,  workpiece  materials  in  the  weld 
region  experience  large  plastic  deformations  during  LFW. 
Under  these  circumstances,  the  use  of  a  pure  Lagrangian 
approach  in  which  the  finite-element  mesh  is  attached  to  and 
moves  with  the  material  may  be  prone  to  serious  numerical 
problems  (due  to  excessive  mesh  distortion).  To  overcome  this 
potential  problem,  an  Arbitrary  Lagrangian-Eulerian  (ALE) 
formulation  is  used  within  which  adaptive  re-meshing  is  carried 
out  of  the  highly  distorted  regions  of  the  workpieces  to 
maintain  a  good  quality  mesh. 

The  fully  coupled  thermo-mechanical  ALE  problem  asso¬ 
ciated  with  LFW  process  modeling  is  solved  numerically  using 
an  explicit  solution  algorithm  implemented  in  ABAQUS/ 
Explicit  (Ref  33),  a  general-purpose  finite-element  solver.  As 
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will  be  shown  later,  a  material  user-subroutine  was  developed 
and  linked  with  this  software  in  order  to  implement  the 
Carpenter  Custom  465  material  model  (presented  in  the  next 
section). 

3.10  Material  Thermo-Mechanical  Models 

For  a  complete  definition  of  the  LFW  process  model, 
specification  of  the  thermo-mechanical  constitutive  models  for 
the  subject  material(s)  is  mandatory.  Unfortunately,  an  over¬ 
view  of  the  public-domain  literature  did  not  reveal  the  existence 
of  such  models  for  Carpenter  Custom  465.  Consequently,  one 
of  the  main  objectives  of  the  present  work  was  the  development 
of  the  thermo-mechanical  material  model  for  this  material. 
Details  of  the  procedure  used  and  the  results  obtained  are 
presented  in  the  next  section. 

3.11  Computational  Accuracy,  Stability  and  Cost 

To  keep  the  computational  cost  reasonable  while  ensuring 
accuracy  and  stability  of  the  (conditionally  stable,  explicit) 
computational  procedure,  a  mass-scaling  algorithm  is  used. 
This  algorithm  adaptively  adjusts  material  density  in  the  critical 
(time-step-controlling)  finite  elements  without  significantly 
affecting  the  computational-analysis  results. 


4.  Material  Model  Development 

In  the  present  analysis,  the  two  workpieces  to  be  welded  are 
assumed  to  both  be  made  of  Carpenter  Custom  465,  HI 000. 
Consequently,  only  one  thermo-mechanical  material  model  had 
to  be  developed  in  the  present  work. 

4.1  Mechanical  Model 

Following  our  prior  work  (Ref  24),  the  mechanical  response 
of  the  workpiece  material  is  assumed  to  be  of  an  elastic/plastic 
character.  The  elastic  response  of  the  material  is  assumed  to  be 


isotropic,  linear,  and  temperature  dependent,  and  to  be  gov¬ 
erned  by  a  generalized  Hooke’s  Law.  As  far  as  the  plastic 
response  of  the  material  is  concerned,  it  is  assumed  to  be  strain 
hardenable,  strain  rate  sensitive,  and  thermally  softenable,  and 
to  be  described  by  the  Johnson-Cook  material  model  (Ref  34). 
Plastic  anisotropy,  as  quantified  by  the  Lankford  coefficient  R, 
is  neglected  since  in  Carpenter  Custom  465,  HI 000,  R  takes  on 
values  near  1.0  (Ref  19). 

The  room-temperature  Young’s  modulus  of  Carpenter 
Custom  465,  HI 000,  has  been  reported  in  Ref  19  as 
199  GPa.  On  the  other  hand,  no  public-domain  data  pertaining 
to  the  effect  of  temperature  on  the  Young’s  modulus  of  this 
material  could  be  identified.  To  overcome  this  lack  of  data,  a 
computational  procedure  is  adopted  within  which  it  is  assumed 
that  temperature  dependence  of  the  Young’s  modulus  in 
Carpenter  Custom  465,  HI 000,  is  identical,  on  a  relative  basis, 
to  that  of  another  related  stainless  steel.  Examination  of  the 
public-domain  literature  identifies  temperature  dependence  of 
the  Young’s  modulus  of  440C  martensitic  stainless  steel  in 


Fig.  6  The  effect  of  contact  pressure  and  surface  roughness  on 
stainless-steel/stainless-steel  interfacial  thermal  conductance  in  air 
(Ref  32) 


Fig.  7  Temperature  dependence  of  (a)  the  Young’s  modulus  and  (b) 
the  Poisson’s  ratio  in  Carpenter  Custom  465,  HI 000,  derived  using  a 
property  correlation  analysis  and  experimental  data  reported  in  Ref  35 


Journal  of  Materials  Engineering  and  Performance 


Volume  23(6)  June  2014—2189 


Ref  35  as  the  best  source  of  data  for  inferring  temperature 
dependence  of  Carpenter  Custom  465,  HI 000.  The  aforemen¬ 
tioned  procedure  yielded  temperature  dependence  of  the 
Young’s  modulus  of  Carpenter  Custom  465,  HI  000,  as  depicted 
in  Fig.  7(a).  Following  the  same  procedure  and  the  same  source 
reference,  temperature  dependence  of  the  Poisson’s  ratio  in 
Carpenter  Custom  465,  HI 000,  has  been  determined,  Fig.  7(b). 

As  mentioned  earlier,  plastic  behavior  of  the  workpiece 
material  during  LFW  is  represented  using  the  Johnson-Cook 
model  (Ref  34).  This  model  is  capable  of  accounting  for  the 
material  behavior  under  large-strain,  high-deformation  rate,  and 
high-temperature  conditions,  the  type  of  conditions  commonly 
encountered  in  the  problem  of  computational  modeling  of  the 
LFW  process.  Within  the  Johnson-Cook  material  model,  the 
workpiece  material  yield  strength  is  assumed  to  be  controlled 
by  strain-  and  strain-rate  hardening,  as  well  as  by  reversible 
thermally  activated  slip-controlled  thermal-softening  effects, 
and  is  given  by  the  following  functional  relation: 

< yy=[A  +  [1  +  C,ln(teyiP‘)]  [l  -  T™} ,  (Eq  3) 

where  ay  is  the  yield  strength;  fpl  is  the  equivalent  plastic 
strain;  fpl  is  the  equivalent  plastic  strain  rate;  €pl  is  a  reference 
equivalent  plastic  strain  rate;  A  is  the  zero-plastic-strain,  unit- 
plastic  strain  rate,  room-temperature  yield  strength;  B  is  the 
strain-hardening  constant;  n  is  the  strain-hardening  exponent, 
Ci  is  the  strain  rate  constant;  m  is  the  thermal-softening  expo¬ 
nent;  and  Th  =  (T-Tvoom)/(TmQlt- Troom)  a  room-temperature 
(C-oom)  based  homologous  temperature  while  Tmelt  is  the  melt¬ 
ing  temperature.  All  temperatures  are  given  in  Kelvin. 

As  stated  earlier,  there  are  no  public-domain  data  pertaining 
to  the  strain-/strain-rate  hardening  and  thermal  softening  of 
Carpenter  Custom  465,  HI 000.  In  other  words,  parameteriza¬ 
tion  of  the  Johnson-Cook  material  model  is  not  available  and 
could  not  be  constructed  based  on  the  available  experimental 
data.  To  overcome  this  lack  of  data,  a  material-property 
correlation  analysis  is  developed  and  applied.  Within  this 
analysis,  it  is  assumed  that  (a)  the  Johnson-Cook  material- 
model  parameters  are  related  to  the  basic  material  properties 
such  as  the  stacking  fault  energy,  the  number  of  independent 
slip  systems,  the  character  of  the  dominant  strengthening 
mechanism  (e.g.,  precipitation-hardening,  strain-hardening,  and 
solution-strengthening),  etc.;  and  (b)  available  Johnson-Cook 
material-model  parameters  for  a  related  material  can  be  used  to 
assess  the  required  parameters  for  Carpenter  Custom  465, 
HI 000.  Specifically,  to  parameterize  the  Johnson-Cook  material 
model  for  Carpenter  Custom  465,  HI 000,  the  following  steps 
were  taken: 

(a)  The  ascending  portion  of  the  room-temperature  uniaxial 
engineering-stress/strain  curve  obtained  under  a  nominal 
plastic  strain  rate  of  0.01  s-1  (Ref  36),  Fig.  8,  is  first  fitted 
to  a  function  oy  =  [A  +  5(ep1)"] ,  using  a  non-linear 
regression  analysis.  This  procedure  yielded  the  following 
Johnson-Cook  material-model  parameters:  A  =  1470  MPa, 
B  =  1575  MPa,  n  =  0.53,  and  t?1 * * *  =  0.01  s-1; 

(b)  The  strain-rate  hardening  parameter,  C,  is  determined 
using  the  following  steps: 


Fig.  8  The  room  temperature  uniaxial  engineering-stress/strain 
curve  for  Carpenter  Custom  465,  HI 000,  obtained  under  a  nominal 
plastic  strain  rate  of  0.01  s-1  (Ref  36) 


r  1  9<*y 

A  d(lntP')  ’ 


(Eq  4) 


(ii)  The  second  term  on  the  right-hand  side  of  Eq  (4)  is 
next  expressed  using  the  conventional  thermally  acti¬ 
vated  deformation  theory  in  which  the  rate  of  plastic 
deformation  is  defined  by  the  following  Arrhenius- 
type  equation: 


tpl 


^P'exp 


(Eq  5) 


where  fpl  is  an  equivalent  plastic  strain-rate  parameter 
(which  scales  with  the  attempt  frequency  of  dislocations 
to  overcome  the  rate-controlling  obstacles  to  slip-based 
deformation),  AG  is  the  associated  activation  energy,  and 
R  is  the  universal  gas  constant.  Within  the  theory  of  ther¬ 
mally  activated  slip,  e.g.,  (Ref  37),  the  activation  energy 
AG  is  defined  by  the  following  expression: 

AG  =  AG0[l  -  (y/f,  (Eq  6) 


where  AG0  is  the  maximum  activation  energy  associated 
with  dislocation  motion  in  the  absence  of  any  applied 
stress,  s  is  the  material  yield  strength  in  the  absence  of  ther¬ 
mal  activation,  and  p ,  q  are  obstacle-shape  defining  param¬ 
eters.  Parameters  AG0,  p,  and  q  are  characteristics  of  (and 
their  values  are  governed  by  the  slip-controlling)  disloca¬ 
tion-obstacle  interactions. 

By  combining  Eq  (5)  and  (6)  and  rearranging,  the  follow¬ 
ing  expression  is  obtained: 


(i)  by  differentiating  yield  stress  as  defined  by  Eq  (3),  in 

which  the  strain-hardening  and  thermal-softening 

terms  are  omitted,  with  respect  to  the  natural  loga¬ 

rithm  of  the  plastic  strain  rate,  it  is  found  that 


After  differentiating  Eq  (7)  to  obtain  doy/d(\mp])  (the 
resulting  expression  is  too  cumbersome  to  be  shown),  the 
resulting  equation  is  evaluated  at  T  =  300  K  and 
epl  =s  0.01  s-1.  During  this  evaluation,  it  is  assumed  that 
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due  to  attendant  high  dislocation  density  in  lath  martensite 
of  Carpenter  Custom  465,  HI 000,  interaction  of  glide  dislo¬ 
cations  with  (repulsive)  forest  dislocations  is  the  deforma¬ 
tion-mechanism-controlling  material  strength.  Conse¬ 
quently,  typical  values  for  this  type  of  dislocation/obstacle 
interactions  are  assigned  to  the  five  obstacle-characterizing 
parameters,  i.e.,  AG0  =  400  kJ/mol,  s  =  2500  MPa, 
epl  =  10  s-1,  p  =  2/3,  and  q  =  2.  The  resulting  value  for 
doy / <9(lnepl)  is  used  in  Eq  (4)  to  obtain  C  =  0.0347. 

(c)  According  to  the  Johnson-Cook  material  model,  thermal 
softening  is  governed  by  two  material  parameters,  TmQ\i 
and  m.  The  Tme]t  =  1673  K  for  Carpenter  Custom  465, 
HI 000,  has  been  reported  in  Ref  19.  Parameter  m  is  not 
generally  assigned  a  physical  meaning  and  is  treated 
mainly  as  a  fitting  parameter.  However,  it  is  generally 
found  that  in  materials  with  comparable  chemistries  and 
micro  structures,  this  parameter  takes  on  similar  values. 
Adopting  this  observation,  m  parameter  for  Carpenter 
Custom  465,  HI 000,  has  been  set  to  a  value  (; m  =  1.75) 
equal  to  an  average  of  values  for  several  precipitation- 
hardened  martensitic  steels  (Ref  36). 

For  convenience,  a  summary  of  the  general,  room-temper¬ 
ature  linear  elasticity,  and  Johnson-Cook  material-model 
parameters  for  Carpenter  Custom  465,  HI 000,  determined  in 
the  present  work,  is  provided  in  Table  1 . 

It  should  be  noted  that  in  the  original  Johnson-Cook  material 
model,  temperature  affects  the  material  strength  only  reversibly 
by  promoting  thermal  activation  of  dislocation  motion.  In  other 
words,  no  permanent  changes  in  the  room-temperature  material 
micro  structure  and  properties  are  assumed  to  result  from  a  high- 
temperature  exposure  of  the  material.  This  assumption  is  not 
fully  justified  in  the  case  of  LFW,  where  it  is  commonly 
observed  that  high  temperatures  in  the  weld  give  rise  to 
dynamic  recovery/recrystallization  with  the  accompanying  sub¬ 
grain  formation  and  grain-size  refinement.  In  addition,  in  age- 
hardened  materials  such  as  Carpenter  Custom  465,  exposure  of 
the  material  to  high  temperatures  may  lead  to  over-aging, 
precipitate  dissolution,  and  austenite  reversion.  Following  the 
approach  proposed  in  Ref  16,  potential  modifications  to  the 
Johnson-Cook  material  model  are  being  investigated  in  our 
ongoing  work.  Specifically,  parameter  A  in  Eq  (3)  is  being 
considered  not  as  a  constant,  but  rather  as  an  instantaneous- 
microstructure-dependent  quantity. 

While  the  original  Johnson-Cook  material  model  is  available 
as  a  built-in  model  within  ABAQU S/Explicit  (Ref  33)  and  the 


Table  1  General,  room-temperature  linear  elasticity,  and 
Johnson-Cook  plasticity  material-model  parameters  for 
Carpenter  Custom  465 


Parameter 

Symbol 

Units 

Value 

Density 

P 

kg/m3 

7840 

Young’s  Modulus 

E 

GPa 

199.0 

Poisson’s  ratio 

V 

N/A 

0.28 

Reference  strength 

A 

MPa 

792.0 

Strain-hardening  parameter 

B 

MPa 

510.0 

Strain-hardening  exponent 

n 

N/A 

0.26 

Strain-rate  coefficient 

C 

N/A 

0.014 

Room  temperature 

T 

1  room 

K 

293 

Melting  temperature 

Tmelt 

K 

1673.0 

Temperature  exponent 

m 

N/A 

1.03 

user  is  only  required  to  provide  values  of  the  parameters  A,  B , 
n ,  C,  €pl,  Tmelt,  Troom,  and  m ,  this  model  is  implemented  in  a 
user-material  subroutine  (VUMAT)  and  linked  with  the  ABA¬ 
QU S/Exp  licit  solver.  This  was  done  in  order  to  initiate 
development  and  implementation  of  the  aforementioned  mod¬ 
ifications  to  the  Johnson-Cook  material  model.  Full  implemen¬ 
tation  of  these  modifications  will  be  presented  in  our  future 
work.  To  validate  the  current  implementation  of  the  Johnson- 
Cook  material  model,  several  LFW  cases  were  analyzed  using 
both  the  user-material  model  and  the  built-in  Johnson-Cook 
material  model  (with  the  same  parameterizations).  The  corre¬ 
sponding  results  (not  shown  for  brevity)  are  found  to  be 
practically  identical,  validating  the  current  implementation  of 
the  Johnson-Cook  material  model. 

4.2  Thermal/Thermo-Elastic  Model 

As  far  as  the  thermal  response  of  Carpenter  Custom  465  is 
concerned,  it  is  assumed  to  be  isotropic  and  temperature 
dependent.  The  thermal  model  used  includes  two  material 
parameters:  (a)  thermal  conductivity,  and  (b)  specific  heat,  cp. 
Temperature  dependencies  of  these  parameters  (Ref  19,  37)  are 
depicted  in  Fig.  9(a)  and  (b).  In  addition,  the  temperature 
dependencies  of  the  true  and  room-temperature-based  mean 
values  of  the  linear  thermal  expansion  coefficient,  a  thermal- 
elastic  parameter  which  controls  the  extent  of  thermal  strains, 
for  the  same  material  are  shown  in  Fig.  9(c)  and  (d)  (Ref  19). 

4.3  Validation  of  the  Material  Model 

To  validate  the  thermo-mechanical  material  model  for  Car¬ 
penter  Custom  465,  HI 000  developed  in  the  present  work, 
simulations  of  a  high-temperature  ring-compression  test  are 
carried  out,  and  the  results  obtained  were  compared  with  their 
experimental  counterparts  (generated  in  our  ongoing  experimen¬ 
tal  investigation).  Within  the  ring-compression  test,  a  ring- 
shaped  test  sample  is  compressed  between  two  flat  platens. 
During  compression,  the  outer  diameter  of  the  specimen 
increases.  On  the  other  hand,  depending  on  the  magnitude  of 
the  friction  coefficient,  the  inner  diameter  may  either  decrease  (in 
the  case  of  a  relatively  large  friction  coefficient)  or  increase  (in  the 
case  of  a  relatively  small  friction  coefficient).  Thus,  the  comparison 
of  the  simulated  and  experimental  ring-compression  results  does 
not  only  help  validate  the  thermo-mechanical  material  model  but 
also  the  contact  model  developed  in  Section  III. 

The  ring-shaped  test  specimen  had  the  following  dimen¬ 
sions:  inner  diameter  =  6  mm,  outer  diameter  =11  mm,  and 
height  =  6  mm.  In  order  to  replicate  the  contact  conditions 
encountered  in  the  LFW  process  simulations,  both  the  test 
specimen  and  the  platens  are  made  of  Carpenter  Custom  465, 
HI 000.  The  results  presented  and  discussed  in  this  section  were 
obtained  at  a  test  temperature  of  973  K,  and  at  a  test-machine 
cross-head  velocity  of  0.2  mm/s. 

In  Fig.  10(a)  and  (b),  a  comparison  is  provided  between  the 
computed  and  experimental  results  pertaining  to  the  evolution 
of  the  inner  and  outer  diameters,  respectively,  of  the  ring- 
shaped  specimen  in  the  course  of  compression.  It  is  seen  that 
the  agreement  between  the  computed  and  experimental  results 
is  fairly  good. 

In  Fig.  10(c),  a  comparison  is  given  between  the  computed 
and  experimental  results  pertaining  to  the  load  vs.  specimen- 
height  reduction  relationship  in  the  course  of  compression. 
Again,  the  agreement  between  the  two  sets  of  results  is 
reasonable. 
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Fig.  9  Variation  of  (a)  thermal  conductivity,  (b)  specific  heat,  (c)  true  linear  thermal  expansion  coefficient,  and  (d)  room-temperature -based 
mean  thermal  expansion  coefficient  with  temperature  in  Carpenter  Custom  465  with  different  heat-treated  conditions  (Ref  19,  34).  Please  see  text 
for  explanation  of  the  line  labels 


In  Fig.  10(d),  a  comparison  is  given  between  the  computed 
and  experimental  results  pertaining  to  the  temperature  evolution 
at  a  point  initially  located  as  follows:  (a)  vertical  position — half 
the  specimen  height,  and  (b)  radial  position — the  average  of  the 
inner  and  outer  diameters.  Again,  the  agreement  between  the 
two  sets  of  results  is  acceptable. 

Based  on  the  findings  reported  in  this  section,  it  appears  that 
the  newly  developed  thermo-mechanical  material  model  for 
Carpenter  Custom  465,  HI 000,  as  well  as  the  contact  model 
provide  reasonable  predictions  of  the  Carpenter  Custom  465, 
HI 000,  mechanical  behavior  under  high-temperature  and  high- 
contact-pressure  loading  conditions. 


5.  Results  and  Discussion 

In  this  section,  a  few  results  pertaining  to  the  application  of 
the  newly  developed  LFW  process  model  to  Carpenter  Custom 


465,  HI 000  (represented  by  the  newly  developed  thermo¬ 
mechanical  material  model),  are  presented  and  discussed.  These 
results  pertain  to  the  spatial  distribution  and  temporal  evolution 
of  the  temperature  and  deformation  fields  (the  fields  which 
affect  micro  structure  evolution  and  spatial  distribution  within 
the  weld  regions).  While  the  prediction  of  the  microstructural 
changes  within  the  weld  is  the  subject  of  our  ongoing 
investigation,  an  example  of  the  results  pertaining  to  the  weld 
micro  structure  will  also  be  presented. 

Specifically,  the  results  presented  and  discussed  in  this 
section  pertain  to  the  effect  of  the  LFW  process  parameters  in 
Carpenter  Custom  465,  HI 000,  on  (a)  spatial  distribution  and 
temporal  evolution  of  the  interface  temperature,  (b)  temporal 
evolution  and  geometry  of  the  expelled  material/flash,  (c) 
temporal  evolution  of  the  axial  shortening,  and  (d)  as-welded 
micro  structure  prediction.  The  LFW  process  parameters  con¬ 
sidered  here  include  (a)  contact-pressure  history  throughout 
different  phases  of  the  LFW  process,  (b)  reciprocation  fre¬ 
quency,  (c)  reciprocation  amplitude,  and  (d)  number  of 
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Fig.  10  Comparison  between  computed  and  experimental  results  pertaining  to  the  evolution,  during  973  K  ring-compression,  of  (a)  inner  diam¬ 
eter,  (b)  outer  diameter,  (c)  load  vs.  specimen-height  reduction  relationship,  and  (d)  temperature  at  a  point  initially  located  at  half  the  specimen 
height  and  at  the  average  of  the  inner  and  outer  diameters 


reciprocation  cycles  or  total  duration  of  reciprocating  motion.  A 
schematic  is  provided  in  Fig.  11(a)  and  (b)  in  order  to  help 
clarify  the  definition  and  time  history  of  the  LFW  process 
parameters.  An  example  of  the  test  matrix  used  in  the  present 
work  showing  the  levels  of  the  reciprocation  amplitude, 
reciprocation  frequency,  and  contact  pressure  in  the  first  three 
LFW  phases  is  displayed  in  Table  2. 

5.1  Spatial  Distribution  and  Temporal  Evolution 
of  the  Interface  Temperature 

From  the  standpoint  of  obtaining  a  sound  weld  with  good 
mechanical  properties,  temperature  uniformity  over  the  contact 
interface,  as  well  as  the  maximum  value  of  the  interface 
temperature,  attained  during  the  course  of  the  LFW  process  is 
of  major  concern.  This  was  the  main  reason  that  temperature 
uniformity  and  maximum  interface  temperature,  during  LFW, 
are  recorded  and  analyzed  in  the  present  work. 


An  example  of  typical  results  pertaining  to  the  temperature 
(in  K)  distribution  over  the  contact  interface  within  the  LFW 
equilibrium  phase  is  displayed  in  Fig.  12(a).  The  results 
displayed  in  this  figure  show  that  the  temperature  difference 
between  the  central  region  of  the  contact  interface  and  the 
interface  edges  is  ca.  35  K.  The  main  reason  for  the  observed 
non-uniformity  in  the  temperature  distribution  along  the  contact 
interface  is  the  periodic  loss  of  contact  between  the  faying 
surfaces  over  edge  regions  which  extend  orthogonally  to  the 
direction  of  reciprocation  and  whose  widths  are  equal  to  the 
reciprocating  amplitude.  In  these  regions  of  the  contact 
interface,  temperatures  are  lower  for  two  main  reasons:  (a) 
heat  is  not  being  generated;  and  (b)  heat  is  being  lost  via 
convection  and  radiation,  while  these  regions  are  being 
exposed,  i.e.,  not  in  contact. 

Depending  on  the  choice  of  the  LFW  process  parameters 
(reported  in  Table  2),  the  maximum  temperature  difference 
over  the  contact  interface  was  found  to  range  between  30  and 
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180  K.  Clearly,  for  improved  quality  of  the  LFW  joint,  LFW 
process  parameters  should  be  selected  in  such  a  way  that  this 
temperature  difference  is  minimized.  The  LFW  process  simu¬ 
lation  parametric  study  carried  out  in  the  present  work  clearly 
demonstrated  that  this  can  be  achieved  through  (a)  reduced 
reciprocation  amplitude,  (b)  increased  reciprocation  frequency, 
and  (c)  increased  levels  of  the  contact  pressure.  This  is  depicted 
graphically  in  Fig.  13,  which  shows,  at  a  constant  level  of  the 
reciprocation  amplitude  (=2  mm),  the  effect  of  contact  pressure 
and  the  reciprocation  frequency  on  the  maximum  temperature 
difference  over  the  contact  interface  during  the  LFW  equilib¬ 
rium  phase.  Due  to  the  aforementioned  difficulties  associated 
with  the  measurements  of  interfacial  temperature  and  its 
temporal  evolution  and  spatial  distribution,  the  results  pre- 
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Fig.  11  Schematic  variations  of  (a)  the  contact  pressure  and  (b)  the 
reciprocating  motion  parameters  throughout  different  LFW  phases 


sented  in  Fig.  12(a)  and  (b)  and  13  cannot  be  validated  directly. 
An  indirect  validation  of  these  results  was  obtained  in  our 
ongoing  experimental  investigation,  which  showed  that  the 
temperature  uniformity  over  the  contact  interface  results  in 
micro  structural  uniformity  over  this  interface. 

The  second  aspect  of  the  interface  temperature  is  its 
maximum  value  attained  during  the  LFW  process.  Clearly,  this 
maximum  temperature  must  not  exceed  the  material’s  solidus 
temperature.  However,  attainment  of  this  condition  may  still 
yield  undesirable  results  regarding  the  mechanical  properties  of 
the  weld  if  the  maximum  interface  temperature  exceeds  the 
temperature  (As)  for  the  onset  of  martensite  — >  austenite 
transformation  (estimated  as  1045  K  for  Carpenter  Custom 
465,  HI 000).  An  example  of  the  typical  results  pertaining  to 
temporal  evolution  of  temperature  at  the  center  of  the  contact 
interface  is  shown  in  Fig.  12(b).  It  is  seen  that,  during  the 
transition  phase,  temperature  steeply  rises  and  that,  during  the 
equilibrium  phase,  temperature  oscillates  over  a  relatively 
narrow  range.  Furthermore,  examination  of  Fig.  12(b)  suggests 
that  there  is  a  possibility  for  formation  of  austenite  at  this 
location  (the  location  which  experiences  the  highest  tempera¬ 
tures  in  the  course  of  LFW),  since  the  maximum  temperature 
(ca.  1100  K)  exceeds  the  As  temperature.  The  aforementioned 
parametric  study  revealed  that  the  maximum  temperature 
experienced  by  the  center  of  the  contact  interface  can  be  made 
lower  than  As  by  proper  selection  of  the  LFW  process 
parameters.  In  addition,  one  must  take  into  account  the  fact 
that  the  material  residing  at  the  contact  interface,  i.e.,  the 
material  that  is  subjected  to  the  highest  temperatures,  is 
continuously  expelled  outward  to  form  flash.  Thus,  formation 
of  austenite  appears  to  be  less  of  a  concern  during  the  LFW 
process.  Furthermore,  due  to  a  continuous  expulsion  of  the 
material  from  the  contact  surfaces,  material  over-aging  within 
the  weld  during  the  LFW  process  appears  to  be  of  less  concern. 
These  findings  have  been  confirmed  in  our  ongoing  weld- 
micro  structure  characterization  experimental  investigation, 
which  confirmed  the  absence  of  melting,  austenite,  and 
excessive  over-aging  in  the  weld  region. 

5.2  Temporal  Evolution  and  Geometry  of  the  Expelled 
Material/Flash 

As  explained  earlier,  LFW  is  accompanied  by  the  formation 
of  flash.  Flash  contains  the  material  which  was  previously 
located  at  the  contact  surfaces  and  was,  in  the  course  of  LFW, 
heated,  softened,  and  subsequently  expelled  in  the  direction  of 
reciprocation.  Formation  of  the  flash  during  LFW  is  highly 
critical  since  the  expelled  material  was,  before  expulsion, 
typically  oxidized,  contaminated,  or  otherwise  compromised 
and  its  removal  creates  clean,  virgin-material  surfaces  with  high 
affinity  for  adhesion/bonding. 

An  example  of  the  typical  results  pertaining  to  the  temporal 
evolution  of  flash  during  the  equilibrium  phase  of  LFW  of 


Table  2  An  example  of  the  Carpenter  Custom  465,  HI 000,  LFW  process  parameter  test  matrix  used  in  the  present  work 


LFW  phase 

Parameter,  units 

Initial 

Transition 

Equilibrium 

Deceleration 

Reciprocation  amplitude,  mm 

0.0 

0.5,  1.0,  1.5,  2.0 

0.0 

Reciprocation  frequency,  Hz 

Contact  pressure,  MPa 

0 

200,  400,  600 

25,  35,  45,  55 

0 

600,  800 
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Fig.  12  Examples  of  the  typical  temperature -based  results  obtained  during  LFW  process  simulation  of  Carpenter  Custom  465,  HI  000,  (a)  typi¬ 
cal  temperature  (in  K)  distribution  over  the  contact  interface  and  (b)  temporal  evolution  of  temperature  at  the  contact-interface  center  point 


Carpenter  Custom  465,  HI 000,  is  displayed  in  Fig.  14(a)  and 
(d).  It  should  be  noted  that  a  surface  temperature  contour  plot  is 
superimposed  onto  the  material  distribution  plot  in  this  figure  in 
order  to  help  relate  the  local  temperature  with  the  extent/rate  of 
flash  formation.  Examination  of  the  results  displayed  in 
Fig.  14(a)  and  (d)  show  that  a  single  flash  (on  each  weldment 
side)  is  formed,  suggesting  a  good  quality  weld.  This  finding 
was  confirmed  in  our  ongoing  experimental  investigation. 


5.3  Temporal  Evolution  of  the  Axial  Shortening 

As  the  thermally  softened  interface  material  is  expelled,  the 
height/thickness  of  the  two  workpieces  being  welded  is 
reduced.  This  phenomenon  is  commonly  referred  to  as  axial 
shortening.  It  is  generally  accepted  that  a  minimum  value  of 
axial  shortening  (which  ensures  complete  removal  of  the 
oxidized,  contaminated  or  otherwise  compromised  material 
from  the  contact  surfaces)  is  necessary  for  the  attainment  of  a 
sound  weld  with  good  mechanical  properties.  It  is  equally  well 
recognized  that  excessive  axial  shortening  leads  to  a  loss  of 


material  and  productivity,  as  well  as  to  an  increased  cost 
associated  with  removal  and  handling  of  excessive  flash. 

An  example  of  typical  results  pertaining  to  the  temporal 
evolution  of  axial  shortening  during  LFW  of  Carpenter  Custom 
465,  HI 000,  is  displayed  in  Fig.  15(a).  For  improved  clarity, 
different  LFW  phases  are  denoted  in  this  figure.  It  is  seen  that 
shortening  begins  in  the  initial  phase  (due  to  application  of  the 
contact  pressure)  and,  after  a  short  transition  phase,  continues 
throughout  the  equilibrium  phase  at  a  fairly  constant  average 
rate.  In  the  forging  phase,  application  of  a  higher  upset  contact 
pressure  causes  an  abrupt  jump  in  the  axial  shortening. 

The  average  slope  of  the  axial-shortening  vs.  time  curve  in 
the  LFW  equilibrium  phase  defines  the  rate  of  axial  shortening, 
also  sometimes  referred  to  as  the  burn-off  rate.  This  quantity 
plays  an  important  role  in  the  LFW  process,  since  it  is  generally 
believed  that  its  value  must  exceed  a  minimum  critical  level  for 
the  attainment  of  a  sound  weld.  The  rate  of  axial  shortening  is  a 
function  of  the  LFW  process  parameters.  This  is  demonstrated 
in  Fig.  15(b)  in  which,  at  a  constant  level  of  reciprocation 
amplitude  of  2  mm,  a  contour  plot  is  displayed  showing  the 
effect  of  the  (friction)  contact  pressure  and  the  reciprocation 
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frequency  on  the  equilibrium-phase  axial- shortening  rate. 
Examination  of  the  results  displayed  in  these  figures  reveals 
that,  as  expected,  both  increases  in  the  friction  contact  pressure 
and  reciprocation  frequency  lead  to  increases  in  the  axial- 
shortening  rate.  The  axial- shortening  and  bum-off  results 
presented  in  this  section  could  not  be  currently  validated.  In 
our  ongoing  experimental  investigation,  efforts  are  underway  to 
develop  axial-shortening/bum-off  rate  measuring  capabilities. 

5.4  As- Welded  Microstructure  Prediction 

As  demonstrated  in  the  previous  section,  the  LFW  process 
model  enables  the  prediction  of  the  effect  of  various  process 


Fig.  13  The  effect  of  the  contact  pressure  and  reciprocation  fre¬ 
quency  on  the  maximum  temperature  difference  (in  K)  along  the 
contact  interface  of  LFW  of  Carpenter  Custom  465,  HI 000,  at  a 
constant  level  of  the  reciprocation  amplitude  of  2  mm 


parameters  on  the  temporal  evolution  and  spatial  distribution  of 
temperature  and  material  flow  fields,  as  well  as  the  time  history 
of  axial-shortening/bum-off  rate.  In  our  ongoing  investigation, 
the  current  work  is  extended  in  two  directions:  (a)  a  module  is 
being  added  to  the  LFW  process  model  in  order  to  enable 
prediction  of  the  material-microstructure  evolution  within  the 
weld,  and  (b)  a  companion  experimental  investigation  of  the 
LFW  process  and  the  weld-micro stmcture  characterization  are 
being  developed.  These  extensions  of  the  present  work  will 
enable  full  validation  of  the  newly  developed  LFW  process 
model.  In  the  remainder  of  this  section,  an  example  of  the 
results,  being  generated  in  our  ongoing  investigation,  which  are 
used  for  the  LFW  process  model  validation  is  presented  and 
discussed. 

In  Fig.  16(a)  and  (b),  a  comparison  is  provided  of  the 
experimental  and  computational  results  pertaining  to  the  spatial 
distribution  of  the  precipitate  mean  radius  (in  microns)  over  the 
mid-plane/contact  surface  of  the  LFW  joint.  The  results 
presented  in  Fig.  16(a)  are  obtained  using  standard  micro- 
structure-characterization  (i.e.,  optical  and  scanning  electron 
microscopy)  techniques.  The  results  presented  in  Fig.  16(b)  are 
obtained  using  the  microstmcture-evolution  module  mentioned 
above.  This  module  enables  determination  of  the  type,  the 
extent,  and  the  temporal  evolution  and  spatial  distribution  of 
micro  structural  changes  within  different  portions  of  the  weld 
[38-41].  The  module  is  based  on  the  key  physical  metallurgy 
concepts  and  principles  of  Carpenter  Custom  465,  HI 000,  and 
includes  the  basic  thermodynamics  and  kinetics  of  various 
interacting  and  competing  phase  transformations  which  may 
take  place  within  the  weld  region  of  this  material.  For  example, 
the  multi-component  nature  of  the  material  is  accounted  for,  as 
well  as  the  kinetics  of  precipitate  nucleation,  growth  and 
coarsening  over  an  arbitrary  local  thermal  history  (the  history 
yielded  by  the  LFW  process  model). 

It  should  be  noted  that  material  residing  within  the  flash  is 
not  accounted  for  in  Fig.  16(a)  and  (b),  for  the  following  two 
reasons:  (a)  this  material  does  not  ultimately  reside  within  the 
weldment  and,  hence,  its  precipitate  micro  structure  is  of  no 


Fig.  14  Typical  results  pertaining  to  the  temporal  evolution  of  flash  during  equilibrium  phase  of  LFW  of  Carpenter  Custom  465,  HI 000, 
obtained  in  the  present  work.  Elapsed  time  between  two  consecutive  field  plots  is  0.12  s.  Please  note  that  a  surface  temperature  contour  plot  is 
superimposed  onto  the  material  distribution  plot  in  order  to  help  relate  the  local  temperature  to  the  extent/rate  of  flash  formation 
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Fig.  15  Typical  results  pertaining  to  (a)  temporal  evolution  of  axial 
shortening  and  (b)  the  effect  of  the  equilibrium-phase  contact  pres¬ 
sure  and  reciprocation  frequency,  at  a  constant  level  of  the  reciproca¬ 
tion  amplitude  of  2  mm,  on  the  rate  of  axial  shortening,  both  for 
Carpenter  Custom  465,  HI  000 

concern;  and  (b)  the  material  residing  within  the  flash  typically 
is  subjected  to  excessive  temperatures  which  may  result  in 
martensite-to-austenite  reversion,  the  phase  transformation 
which  was  not  accounted  for  in  the  microstructure-evolution 
module. 

A  comparison  of  the  results  displayed  in  Fig.  16(a)  and  (b) 
reveals  that  (a)  the  precipitate  size  is  generally  non-uniform 
over  the  contact  interface.  Specifically,  the  largest  precipitate 
sizes  are  found  in  the  innermost  portion  of  the  corresponding 
x-y  section,  the  portion  which  experiences  the  highest  temper¬ 
ature.  The  edge  portions  of  this  section  contain  smaller  particle 
sizes  due  to  the  operation  of  heat-transfer  processes  between  the 
workpieces  being  welded  and  the  surroundings,  (b)  the 
experimental  and  computed  results  pertaining  to  the  spatial 
distribution  of  the  precipitate  mean  radius  over  the  faying 
surfaces  are  in  decent  agreement,  and  (c)  the  overall  size  of  the 
precipitates  is  fairly  well  accounted  for  by  the  micro  structure- 
evolution  module. 


Fig.  16  Comparison  of  (a)  experimental  and  (b)  computational 
results  pertaining  to  the  spatial  distribution  of  the  precipitate  mean 
radius  (in  microns)  over  the  mid-plane/contact  surface  of  the  LFW 
joint  in  Carpenter  Custom  465,  HI 000 


6.  Summary  and  Conclusions 

Based  on  the  work  presented  and  discussed  in  the  present 
paper,  the  following  main  summary  remarks  and  conclusions 
can  be  made: 

1.  A  LFW  process  model  for  Carpenter  Custom  465,  HI 000 
(a  precipitation-hardened  martensitic  stainless  steel),  is 
developed  using  an  Arbitrary  Lagrangian-Eulerian  finite-ele¬ 
ment  framework.  The  model  accounts  for  all  the  key 
mechanical  and  thermal  aspects  of  the  process  and,  as  such, 
is  suitable  for  use  in  a  process-optimization  analysis. 

2.  Since  there  is  no  publicly  available  material  model  for 
Carpenter  Custom  465,  HI 000,  a  thermo-mechanical  con¬ 
stitutive  model  for  this  material  had  to  be  developed, 
parameterized,  and  validated.  This  was  accomplished  by 
(a)  using  all  publicly  available  thermo-mechanical  data 
for  this  material,  (b)  employing  a  material-property  corre¬ 
lation  procedure,  and  (c)  comparing  the  computed  results 
for  a  high-temperature  ring-compression  test  with  their 
experimental  counterparts.  This  approach  confirmed  that 
the  newly  developed  thermo-mechanical  model  for  Car¬ 
penter  Custom  465,  HI 000,  can  reasonably  account  for 
the  mechanical  response  of  this  material  under  high- 
temperature  and  high-contact-pressure  conditions,  the 
type  of  conditions  encountered  in  the  LFW  process. 

3.  Application  of  the  LFW  process  model  clearly  revealed 
the  effect  of  the  key  process  parameters  such  as  reciproca¬ 
tion  frequency,  reciprocation  amplitude,  and  the  contact 
pressure  on  the  spatial  distribution  and  temporal  evolution 
of  the  temperature  and  the  flow  fields  within  the  weld. 
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4.  The  resulting  thermal  and  flow  fields  are  used  within  a 
weld-micro  structure  module  to  predict  the  spatial  distri¬ 
bution  of  the  precipitate  size  over  the  contact  surface. 
The  predicted  results  are  found  to  be  in  decent  agreement 
with  their  experimental  counterparts. 
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